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Self-similar Euler singularities may be useful for understanding some aspects of Navier-Stokes turbu- 
lence. Here, a causal explanation for intermittency is given, based on the control of the sudden growth 
of the gradients by the Euler equations. This explanation uses certain Euler solutions as intermediate 
asymptotics in Navier-Stokes turbulence [Q - controlling the dynamics over a limited spatial and temporal 
^vq domain. These arise from an analysis of similarity equations, previously discussed by Pelz and Green 

which yield experimentally testable predictions. Three main points are presented here: scalings of suitable 
characteristic lengths with time from a critical time I ~ (t a — t) a , a > 1, a discussion of invariant sets of 
the similarity equations that result, and a discussion of cutoff mechanisms. The value a = 3/2 appears to 
correspond to Kolmogorov scaling for turbulence. Some limited experimental evidence is presented from 
Euler ian gradient measurements at the Kolmogorov scale showing 1 < a < 3 values. Much testing is nec- 
essary to ascertain the final usefulness and validity of these ideas, as several conceptual obstacles remain. 
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1 Introduction 

> ' 

[- — , ■ Turbulent flows of liquids and gases occur in a host of natural and engineering phenomena. A detailed 
understanding of atmospheric flows on Earth and the planets, predictive design for industrial equipment, 
and understanding of the processes within the liquid cores of planets all rely on an understanding of 
: turbulence. One might naively think this an easy task, given that the governing equations are well known 
(i.e., the family of equations stemming from the Navier-Stokes equations with appropriate body forces). 
That understanding is hampered by the inherent nature of turbulence itself, where the velocity field is rough 
and velocity gradients are intermittent. Phenomenological considerations of Richardson and Kolmogorov 
have been very fruitful [3J E] , but have also left deep holes in our understanding of the intermittent nature 
of turbulence By intermittency, I mean a host of interlocking observations that the basic observables 
in turbulent flow; the vorticity, dissipation, helicity, acceleration, and their ability to effect advected scalar 
and vector fields (concentration or perhaps magnetic fields) all show abnormal statistics. As one example, 
in experiments one can observe accelerations hundreds times the means, lying 30 standard deviations from 
.£h those means jHEl) or observe the dissipation locally becoming one hundred times its average value [8|l9llT0]. 
^ ■ The probability distributions for these observable then show far from Gaussian statistics with long tails. 
These types of gradient or time derivative concentrations appear local in space and time. Recent results of 
Gibbon and Doering suggest that the time intervals where singularities are possible become more localized 
as the Reynolds number increases 

There is little in the phenomenological picture of a cascade of energy to smaller scale (envisioned in 
wavenumber space) to explain intense events. The framework discussed here is an attempt to explain these 
observations, although it is as yet a hypothetical explanation. Rather than attempting a phenomenological 
explanation, instead I seek a causal explanation rooted in the known equations of motion and the relevant 
forces. 

This is a non-statistical approach, but with the need to ultimately explain the observed statistics. First, 
I argue that the viscous forces are not responsible for the sudden growth locally of gradients, and that in 
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this context, during the growth, one should consider the viscous forces irrelevant. This limits our scope 
to the growth phase of these events, and lands us in the context of the Euler equations. Solutions to the 
Euler equations may explain the sudden growth in the gradients and time derivatives in the context of an 
intermediate asymptotic j^. These solutions will be locally relevant for only a range in times, nearby the 
sudden collapse, but limited to times before viscosity, or perhaps cavitation, act to limit the strength of 
these near-singularities. As an example, one can look to the Burgers equation ^2]. There, negative slopes 
tend to collapse to form shocks, but the gradients are limited by the action of viscosity. In the case when 
the viscosity is quite small, there is a substantial range in time where the gradient value at the location of 
minimum gradient is approximately governed by the inviscid Burgers equations, with the gradients growing 
as predicted, i.e. inverse in time from a critical time. 

What is unlikely to be understood through this exercise are the important long lived vortices so many 
have remarked on which form an independent, but perhaps causally connected part of turbulent 
phenomena. Remarks on the connection between long lived vortex filaments and the local collapse of 
gradients studied here are given in the conclusion. 

When entering the study of gradient focusing and local intermittent events in turbulence, one treads 
on the long worn and active path of singularities in Navier-Stokes and Euler equations |27| I28| I15[ I29| \',VA\ 
EU EH1 Oil EH EH| ■ Although these studies are conceptually related, I view the Euler solutions discussed 
here (which are singularities of the Euler equations) as merely intermediate asymptotics for Navier-Stokes, 
as remarked earlier, and not observable singularities of viscous fluid flow. Many previous studies have 
concentrated on the Leray scaling a = 1/2 |13| |2j or on the case where a co-collapsing sphere shows 
energy conservation a = 2/5 [Tl|- At least one simulation of R.M. Kerr has shown some evidence for 
a = 1 behavior |15| . I restrict my attention to the case (described later) of a > 1 so that, as observed 
in turbulence, the velocities are not growing abnormally large. In many turbulent flows the probability 
distribution of the velocity is close to Gaussian, while the gradients are intermittent. 

Whether or not there are singularities in Navier-Stokes solutions is perhaps of little practical conse- 
quence - extra phenomena such as dilations, cavitation or extreme rarity render them of little relevance. 
One should not mistake the irrelevance of singularities of the Navier-Stokes equations for an argument for 
the same for Euler singularities. Here, I'm proceeding based on the premise that the causal structure of 
Euler singularities is the underlying cause of intermittency in Navier-Stokes turbulence, even if they are 
capped by viscosity or other extra constraints. 

While the following calculations and observations are presented with the hope of better understanding 
Navier-Stokes turbulence, their final utility is yet unknown. Nevertheless, because of the possibility of their 
importance, it seems useful to present them. 

2 Similarity methods 

There have been significant successes in understanding singularities and self-focusing in strongly nonlinear 
physical systems using similarity methods pJElEHEIlEn]- Our previous work on the collapse of waves and 
the production of jets on a free surface illustrates the utility of similarity techniques in allowing analytical 
and numerical progress on systems not otherwise amenable to analysis |21j . While there is no guarantee of 
the success of this approach for Euler collapse, proceeding with this hope may lead to new insight. Observe 
that in local gradient collapse the length scales become small as one approaches maximum focusing. In 
what follows I will assume that useful information can be obtained by consideration of the hypotheses that 
the length scales in all directions behave equally, and that the functional form is a power law in time. 
Regarding the first hypothesis, one might say that the different directions may have trouble behaving in 
different manners as they are constantly effecting each other and being mixed by motions and rotations. 
This argument, however, is hardly compelling. One might also argue that the relevant equations balance 
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linear and nonlinear terms, so power-law time dependences may arise. This, however, is not necessary 
and depends on the algebraic form of these balances. There might be collapse events of more complicated 
forms; understanding their analytical structure may be much more difficult. 
Beginning with the Euler equations for incompressible flow: 

d t v+ (y ■ V) v + ivP = (1) 

V-u = 0, (2) 
one can apply the above assumptions to form an intelligent guess for the form of the solutions: 

if} t \ ( tjc ( t 



P (l -, i) = p( „. )2 (^) M n(A(J_)") w 

Here various quantities are made dimensionless using some characteristic large scale L* , velocity scale v*, 
and time scale t* = L* /v* , as would be natural in the context of a shear flow. One can then seek solutions 
for the time independent similarity functions G and II. The above ansatz assumes all lengths scale as 
I ~ (to — t) a , velocities as v ~ (t — t) a ~ l . It is important to require that a > 1 in order that the velocities 
remain bounded as we approach t Q from below. By defining a scaled coordinate r and a scaled gradient 



x ( r 



L* \t -t 



(5) 



V=(^)Vv (6) 

one can eliminate time dependence from the Euler equations and obtain what are referred to as similarity 
equations 

(l-a)G + a(f-V)G+(G-V)G + Vn = 0, (7) 

V • G = 0. (8) 

Solutions to Eqs. [7| and |H] correspond to solutions to the original Euler equations. Symmetries of the 
Euler equations may be used to generate further solutions associated with the similarity ansatz. The time 
or spatial origin of the solutions may be shifted, and one may also add a constant velocity. 

The solutions to Eq. [7| are assumed to be smooth - they represent solutions which retain their smooth- 
ness up to the critical time t Q , and loose it by the length scales limiting to zero. At that critical time the 
solution remains smooth at all locations (governed by the far field of Eq. EJ) except at the origin, where the 
gradients become unbounded. Pelz and Green [2j have examined the stability properties of these solutions. 
Here I present new results and discussion on the form of the solutions of Eq. their correspondence to 
Lagrangian dynamics, viscous cutoff mechanisms, and their relationship to turbulence with Kolmogorov 
statistics. 

I examine the similarity equations progressively in one, two and three dimensions, in order to better 
understand their function and utility. 
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3 One dimensional similarity equations 



In one dimension, one neglects the requirement of incompressibility and drops the pressure term in Eq. [7J 
In this case: 

(1 — a)G + axd x G + Gd x G = 0. (9) 
This equation has several families of implicit analytical solutions [22] , only one of which is relevant here: 

x = -G + aG a/( - a ~ 1] (10) 

The far field (i.e., large \x\) limit of Eq. ^TOJ is G = x 1 ^ 1 ^. This similarity solution gives the generic 
formation of the first shock in the inviscid Burgers equation, and accurately describes the formation of a 
near-shock in the diffusive case, when the diffusion is small. In particular, Burgers equation is unstable 
where the slope has a maximum negative value, the magnitude of the slope following a characteristic 
diverges as (t a — t)^ 1 for a time until the diffusive cutoff takes over. The value a = 3/2 is generically seen, 
due to a curious reason, and leading to a G ~ \x — a^ol 1 ^ 3 far field. The location for the characteristic which 
sees the first shock in the inviscid case is the most negative slope in the initial condition. Generically, 
the second derivative will vanish there; locally the curve will be dominated by the linear and cubic terms. 
These cubic terms yield, in the confluence of characteristics in the shock, a functional form at the time 
of the shock u = u Q + v* \ (x — x a ) / L*\ 1 /^ . The probability distribution of negative gradients in randomly 
forced Burgers turbulence is governed by these pre-shocks [23] . 



4 Two dimensional similarity equations 

The situation in two dimensions is quite different. Due to the simple advection of vorticity in two dimensions 
|24j - the vorticity remains bounded and no similarity solutions exist. However, examining this in more 
detail will assist in understanding the three-dimensional case. 

In two dimensions, the vorticity vector lies normal to the plain of motion. Due to this, there are no 
gradients of the vorticity along the direction of the vorticity, and one can write the vorticity similarity 
equation (curl of Eq. EJ) as: 

u + a{r ■ V)lo + {G- V)cj = (11) 

I can demonstrate that these equations have no solutions using several observations. First, define the 
associated vector field H = ar + G. This field is used to rewrite the vorticity similarity equations in two 
dimensions as: 



u + (H ■ V)uj = (12) 
It is useful to define the dynamical system: 

which yields characteristics for the vorticity. In terms of the trajectory f(s) one can write the vorticity 
equation along these characteristics as 

d7 = (14) 

The dynamical system for f(s) is everywhere expanding, as V • H = 2a and a is positive. Nearly all tra- 
jectories starting out near the origin escape to infinity. Due to the (assumed) smoothness of the dynamical 
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system, Eq. 1131 there must be an invariant set near the origin which does not escape. As this calculation is 
restricted to two dimensions, the possible invariant sets (by the Poincare-Bendixon theorem |25j ) may be 
fixed points or closed orbits. Along any characteristic r(s) the vorticity must be u>(s) = uj(0)e~ s . At fixed 
points, by definition, H = 0, giving uj = from Eq. ^]at those points. For closed orbits, the vorticity uj 
must also be zero so that after going around one orbit it returns to its original value. Points starting nearby 
the closed orbits or fixed points escape to large radius. As their initial vorticity must be zero, as they move 
out from near the origin their vorticity remains zero. From this one can understand that the vorticity is 
zero everywhere in the plain. The far field (large \x\) solution for G must have a dependence rising like 
r ( a ~ 1 )/ a ; where the first two terms of Eq. 0are dominant. As the vorticity is zero (i.e., V x G = 0), the 
vector field G = V0 must be a gradient field with V 2 </> = 0. The far field dependence is incompatible 
with a smooth Harmonic function (j) which must have an integer n far field dependence r n . The far field 
exponent (a — l)/a is not an integer for all cases of interest a > 1. Thus, the only solution to the similarity 
equations in two dimensions are G = uj = 0. Another proof is also possible using integral methods |26| . 

5 Three dimensional similarity 

equations - Lagrangian dynamics 

The flow generated by H is also of considerable use in understanding the three dimensional case. Trajec- 
tories of Eq. are directly associated with particle trajectories x(t) in lab coordinates. The invariant 
sets are especially important and correspond to trajectories which undergo the most singular Lagrangian 
motions. Trajectories in the laboratory frame follow a path dx/dt = v; this equation is equivalent to Eq. 
blunder the identifications of s = — ln[(i — t)/t*] and Eq. 5 for r. The correspondence between these 
different quantities is summarized in the diagram: 

t — > X 

-ln[(t„-t)/t*] | i x[t*/{t -t)] a /L* (15) 

H 

s — > r 

Note as (t — i) — ► 0, that s — > oo. As the flow generated by H has divergence everywhere positive 
(V ■ H = 3a), trajectories all flow outward toward large radius, and there must be repelling invariant sets 
in the vicinity of the origin (assuming of course that the similarity solution is smooth and exists!). Several 
types of invariant sets are possible: fixed points, closed orbits, or perhaps strange saddles 1 . Which type 
of invariant set will be determined by the precise form of the outer solution. Assuming these solutions 
exist, and pending a numerical solution, it is useful to explore the consequences of these invariant sets. 
Fixed points of the dynamical system Eq. ^Jare defined by H = 0; we label them r*. These points have a 
number of properties and yield testable prediction for the Lagrangian motion. At points f*, one has that 
af= —G = — VII from the definition of H and the similarity equation (Eq. EJ) with H = 0. 

We can see why the non-existence proof for two dimensions fails in three by examining the vorticity 
equation in three dimensions (the curl of the similarity Eq. |7J| : 

uj + (H • V)u = (oj ■ V)G. (16) 

Rather than the vorticity being zero at the fixed points of Eq. I13| the vorticity must be a eigenvector of 
the gradient matrix with unit eigenvalue M6j = uj. Although the 2-D non-existence does not carry over to 
3-D, I have yet been unable to show existence directly, numerics may give some evidence in future work. 

1 2-tori are excluded, since their interior volumes would be invariant, contradicting V ■ H = 3a > which implies that all 
volumes grow exponentially as exp(3as) following the flow defined by Eq. If 31 
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From the equivolency of flows of H and Lagrangian trajectories, one can calculate the form for ob- 
servables due to the fixed points and limit cycles of the similarity flow. For a fixed point r* transformed 
to laboratory coordinates the trajectory is x = L*[(t — t)/t*] a r* + v Q t + x* , where we have allowed for a 
background velocity and origin. The Lagrangian velocity would be x = —aL*[{t — t)/t*] a_1 f*/t* + v ; 
acceleration would be x = a(a — l)L*(t a — t) a ~ 2 r* / it*) 2 . This acceleration diverges at the critical time, 
clearly in need of a cutoff due to viscosity or other physical intervention. These forms, as intermediate 
asymptotics are testable predictions for Lagrangian experiments and numerics - are events observed which 
hold such forms? 

In the case of an invariant cycle in the flow of H one obtains a different type of motion. This rotating 
motion in similarity coordinates, with some characteristic frequency in s, yields an accelerating chirp of 
spiraling motion. For simplicity, let's say some some similarity coordinate has harmonic motion r% = 
Acos(£ls). Carried to laboratory coordinates the motion is 

x 1 = i4L*(^)«co8[nin(^)]. (17) 

The laboratory Lagrangian frequency thus rises as (t a — i) -1 while the radius collapses to zero as (t a — t) a . 
Such motions may be related to the large accelerations observed in recent Lagrangian tracking experiments 
013) and may account partially for the extreme intermittency seen in the acceleration distributions. It is a 
likelihood that a clear understanding of the viscous cutoff is needed to fully account for such distributions; 
some aspects of the viscous cutoff are discussed in Sec. 8. These calculations give a plausible specific cause 
for extreme Lagrangian accelerations - which is independent of the intense long-lived vortices which form 
an alternative explanation for those observations. 

6 Far field behavior 

The large r behavior is particularly significant, as close to the critical time, the near field character has 
shrunk (in laboratory coordinates) into an insignificant volume, leaving a relatively large area displaying 
the far field behavior in its wake. For large r, one can calculate the behavior of three dimensional solutions 
of Eq. [Tj Only the first and second terms are relevant, the nonlinear term and pressure are sub-dominant. 
In that case the resulting equation (1 — a)G + a{r ■ V)G = is analytically solvable: 

G = A a - 1 ^ a f(e,ct>) (18) 

for spherical coordinates r, 8, <p. The function / is a general S 2 — > 1Z 3 vector field subject only to V-G = 0. 

Thus observe that the far field radial dependence is fixed by a, while the angular dependence is general. 
It has been conjectured ,30 that Eq. 0has solutions for all admitable far field forms, while it is important 
to emphasize solutions (numerical or analytical) have not yet been found. 

The far field behavior has particular significance if one hopes to use these solutions to understand 
turbulence. The far field solution when a = 3/2, if applied in many events, would lead to velocity difference 
statistics associated with Kolmogorov scaling and the known scaling for the third order structure function. 
When a = 3/2 the far field behaves as ||G|| ~ r 1 / 3 ) which is retained in the lab frame velocity near 
the critical time. This connection in particular has been a motivation for seeking a connection between 
these Euler similarity solutions and Navier-Stokes turbulence. Why the a = 3/2 case may be prevalent is 
precisely unknown, but discussed in the conclusion. 

7 Near field behavior 

One can learn about the structure of the gradients at the fixed points by examining the similarity equation. 
Taking the gradient of that equation, to obtain a matrix equation for M = diGj, one obtains at the fixed 
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points H = that: 



M + M + n 2 = 0. 



(19) 



This matrix Riccati equation is related to the quadratic nonlinearities in the Euler equations. Gibbon 
[glj has shown how quaternions can be used in this context to relate Euler equations and the equations 
of ideal magneto-hydrodynamics to a Schrodinger equation. Equation 20 yields nine constraints on the 
components of M and H2. Rotating M into a form where the symmetric part is diagonal, and given that 
II2 is symmetric, there are eleven independent components between them. Solving Eq. yields three 
types of solutions, each in two parameters. Further solutions come from permuting the axis in this form, 
or by rotations relaxing the diagonalization of the symmetric part of M. 
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The gradient matrix at these fixed points gives rise to rapidly growing gradients in laboratory coordi- 
nates. At these points, taking the gradient of Eq. 3, yields: 



dxj 



Mv* 



t -t 

t* 



(23) 



leading to the vorticity and strains all growing jointly and suddenly. 



8 Viscous cutoff 

One can obtain estimates for the small scale cutoff for these events by analyzing the relevant forces. For 
all events with a > 1 viscosity will be relevant at some late time, where the length scales become small 
enough for the local Reynolds number to be of order unity. Comparing the inertial forces with viscous 
forces for the scaling ansatz, (u • V)n < uV 2 u yields an expression for the time when (for smaller (t — t)) 
viscosity will be dominant. This happens when 

^ <(ir) i/(i- 2Q)) (24) 
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Figure 1: Time trace of the dissipation measured at the Kolmogorov scale in an oscillating grid experiment. 

where R* = u*L* jv. This leads to a definition of a crossover time tk where the equality holds above. At 
the crossover time, the length scale is 

l k = L*(ir) a/(1 - 2a) ; (25) 

this corresponds to the Kolmogorov scale for that event, and depends on the value of a. In particular, 
events with a = 3/2 have the normally defined Kolmogorov scale, while events with a < 3/2 have a reduced 
minimum scale. This expression is similar to those obtained for the Kolmogorov length in the multifractal 
description of turbulence, where ~ for the Holder exponent h |41j . 

This fluctuating Kolmogorov scale effects the maximum gradient observed in different events. Assuming 
that the crossover to viscous dominating times signals the end of the growth of the gradients, one can 
estimate that 

(yi^)k-^{R*) l,{2a ~ 1 \ (26) 

which depends on a. This estimate serves for both strains and vorticity. The crossover acceleration then 
scales as the maximum value that the inertial term takes 

afc ^ (!f!l!(jj*)(«-2)/(2«-i). (27) 
L* 

This leads to a prediction that the index a effects the maximum acceleration in a particular way. In the 
case a < 2, the crossover acceleration shrinks with increasing Reynolds number, while for a > 2, the 
crossover acceleration grows with acceleration. Presumably the viscosity caps the maximum acceleration, 
but these values still grow without bound with increasing Reynolds number. Note for all of these estimates 
a = 1/2 corresponding to the Leray solution |13| yields an uncapped singularity. 

Another way that the maximum gradients may be capped comes about due to the outer solution. As 
these events must live in a finite domain, and must be patched onto realistic boundary conditions, the 
ideal form cannot exist out to large radius. Thus either through this departure in the far field, or by the 
outer condition deviating from the ideal form, the inner solution may miss the maximally focused condition 
before viscosity would cap the event. 

9 Experimental observations 

Although the main parts of this paper are analytical, with some turbulent gradient data in hand, I present 
some limited support for these ideas. A novel optical instrument has been utilized to measure the gradients 
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Figure 2: The rise in the dissipation for the event shown in Fig. 1 has a power law form for some range 
in times. Many such time series are examined to obtain events when the rise is of power law form for at 
least one decade in time, with a coefficient of regression for this decade of at least R = 0.98. 

diUj(t) at the Kolmogorov scale (see jH] for details about the flow and instrument). Focusing on the most 
extreme gradient events, I presume that one needs to be both looking at the center of the event and near 
the critical time in order to have extremely large gradients. Since the measurement is Eulerian, solutions 
to Eq. 7 will both develop in place, and more probably, be swept by the measurement volume. As these 
are swept by, the time dependence reflects the spatial dependence of the field; near the critical time this is 
the far field dependence. 

The far field has velocities scaling as Vi ~ r^ a_1 ^ Q so that gradients behave as di,uj ~ j-- 1 /". The 
dissipation per unit mass e = (v/2) J2ij{di u j + dj u i) 2 ~ r~ 2 / Q . Fig. 1 shows one example of a time trace 
for a large dissipation event. As events are swept by one might occasion to see e(t) ~ t~ 2 l a . An extensive 
data set of gradients has been examined for events where the dissipation behaves as 

e(t) ~ (t - t) s (28) 

for some t\ < t < ti < t Q . So that they are at least minimally serious scaling I require at least one decade 
of scaling, i.e. (t± — t )/(t2 — t Q ) = 10 with a coefficient of regression for a fit to Eq. EHl R 2 > 0.98. The 
event shown in Fig. 1 is redrawn around the large dissipation growth showing (Fig. 2) one example of 
such a scaling. Figure 3 shows a histogram of the a = —2/5 values for all 892 events observed of this type. 
Note that the Kolmogorov value a = 3/2 occurs near the observed maxima. These observations are only 
intended to partially support the analysis discussed as the bulk of this paper. Further experimental testing 
will be discussed elsewhere so as to not distract from the main points presented here. 
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Figure 3: Histogram of the observed values of the exponent a, measured from the rise in dissipation, 
conditioned as indicated in Fig. 2. 

10 Conclusions 

The observation of events with 1 <= a < 3 may be associated with the much prior commented on 
multiscaling of the velocity field ^B]- The solutions corresponding to K41 statistics is a = 3/2 Why 
these solutions are dominant may be of similar reasons to the Burgers case, where the same exponent 
also dominates. One might hypothesize that locations where the gradients obey M + M 2 = H2 and the 
second derivatives (didjVk) vanish may be the locations susceptible to this types of blowup - in analogy 
to inviscid Burgers shock development. The dominance then of the cubic terms in the velocity near field 
leads, through a characteristic collapse, to a = 3/2 events. While other types of a events also appear 
possible, a selection mechanism is not obvious. 

And then one must deal with the fact that these solutions are self-similar around one point in time and 
space, while turbulence is a statistically self-affine (or multi-affine) phenomena [HI El El- What is needed 
is a number of events associated with Eq. occurring at many spatial locations and times. This may 
be accommodated if these solutions give rise to the growth of additional events nearby due to the form 
the intermediate field takes. In particular, if a sufficiently complicated angular structure for the far field 
occurs, when that must flow into the relatively simple near field, there may be locations with the right 
type of structure (perhaps in the second order derivative field) which gives rise to new growing collapse 
events. 

Finally, these events are distinct from long lived vortices others have observed, but perhaps causally 
connected. A sudden growth of the gradients around a point might be expected when two tube-like vortices 
collide obliquely |32j. Also, I acknowledge that more complicated solutions, perhaps with less trivial length 
scaling may also play an important role. 
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